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Abstract 

In this paper we present a complete linear synchrobetatron coupling formalism by studying the 
transfer matrix which describes linear horizontal and longitudinal motions. With the technique 
established in the linear horizontal-vertical coupling study [D. Sagan and D. Rubin, Phys. Rev. 
ST Accel. Beams 2, 074001 (1999)], we found a transformation to block diagonalize the transfer 
matrix and decouple the betatron motion and the synchrotron motion. By separating the usual 
dispersion term from the horizontal coordinate first, we were able to obtain analytic expressions of 
the transformation under reasonable approximations. We also obtained the perturbations to the 
betatron tune and the Courant-Snyder functions. The closed orbit changes due to finite energy 
gains at rf cavities and radiation energy losses were studied by the 5x5 extended transfer matrix 
with the fifth column describing kicks in the 4-dimension phase space. 
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I. INTRODUCTION 



The synchrobetatron coupling (SBC) comes from dispersion at rf cavities and the path 
length dependence on the amplitude of betatron motion. The dispersion at an rf cavity 
makes the longitudinal kicks received from the cavity affect the betatron motion. Since 
the longitudinal kicks depend on the arrival time of the particles, the longitudinal motion 
is coupled to the betatron motion. On the other hand, particles with different betatron 
amplitudes have different path lengths which affect the arrival time. So betatron motion is 
also coupled to the longitudinal motion. 

Traditionally many authors treated SBC with the Hamiltonian dynamics approach 

HQ, 

which is a general and complete description and naturally covers effects of nonlinearities. It 
is very useful for the study of synchrobetatron resonances since in such cases one can focus 
on only the resonant term of the synchrobetatron potential. However, the Hamiltonian 
approach is cumbersome for the off-resonance cases which are most common for storage ring 
operations. In the linear case, a parallel approach is the matrix formalism first proposed 
by Chao in Ref. 3| which described the construction of the 6x6 transfer matrices and the 
decomposition of the coupled motion to the eigen-modes of the one-turn transfer matrix. 
Ref. p] also described an iterative procedure to include the nonlinear effects. Chao later 
applied the method to compute the beam tilt angle between x-z plane due to rf cavities {4]. 

Recently the study of low-alpha lattices stimulated Shoji's work on the path length effect 
which yielded an important result of bunch lengthening due to betatron emittance and 

n , , 

dispersion [5j. On the other hand Ref. |6J studied the SBC- induced closed orbit change by 
considering the dipole-like kicks in the horizontal betatron phase space due to the sudden 
changes of energy at a nonzero-dispersion rf cavity. The authors derived the horizontal 
closed orbit changes induced by the finite energy gains at the rf cavities and verified with 
both simulations and experiments. 

In this paper, we will study the linear synchrobetatron coupling under the transfer matrix 
framework without considering diffusion and damping due to radiation. Since the vertical 
motion is not essential to the SBC, we don't consider it for simplicity reasons. We then 
study the 4x4 horizontal-longitudinal transfer matrix in the same manner as the horizontal- 
vertical coupling is studied [7|. Namely, we try to decouple the horizontal and longitudinal 
motions by using a coordinate transformation to block diagonalize the transfer matrix and 
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obtain the normal modes, in this case, the pure betatron mode and the pure synchrotron 
mode. We first study the fixed-energy case in which no rf cavity exists (or the rf gap voltage 
is set to zero). The transformed coordinates include the usual betatron coordinates and 
momentum deviation coordinate. But the longitudinal phase coordinate is modified by a 
term involving Dx'p — D'xp which corresponds to the bunch lengthening effect studied by 
Shoji jlj]. In cases with rf cavities, we first apply the previous fixed energy transformation 
to separate the dispersion term. Since the synchrotron motion is usually slow, the coupling 
(off-diagonal) blocks of the transfer matrix for the new coordinates are small. Therefore 

n 

we can perform the block diagonalization procedure proposed in Ref. [7] approximately yet 
with high precision. The transformation matrix is expressed analytically with the usual 
parameters such as the Cour ant- Snyder parameters, dispersion functions and the rf voltage 
slope. When the normal modes are obtained, we can calculate their contributions to the 
beam width and bunch length. Perturbations to horizontal betatron motion due to SBC, 
including changes to the betatron tune and Courant-Snyder functions are also obtained. 
It is well known that SBC causes changes to the beam orbit jf], [si, 1(], H|. In this study 



we find the closed orbit in the 4-dimension phase space in an analytical form through the 
5x5 extended transfer matrix method |3j, with the fifth column containing the 4-dimension 
kicks the beam receives from rf cavities and dipole magnets. Both finite energy gains at rf 
cavities and radiation energy losses in bending dipoles are considered. The radiation energy 
loss is random by nature. However, since usually hundreds of photons are emitted in one 
revolution, much more than the number of dipole magnets in which emission happens, we 
consider the radiation energy loss as a steady and uniform process. The radiation energy 
losses contribute additional terms to the horizontal closed orbit change. The above results 
are verified with the accelerator modeling code AT 

This paper is organized as follows. Section Q] is this introduction. Section [Til describes 
the block diagonalization of the 4x4 transfer matrix in the fixed energy case. Section [TTT1 
presents the matrix formalism of the synchrobetatron coupling. Section HVl is the calculation 
of closed orbit changes induced by finite energy gains at rf cavities and radiation energy 
losses. Section [V] shows simulation results and the comparison to the theory. Section IVII 
gives the conclusions. 
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II. BLOCK DIAGONALIZATION FOR A FIXED ENERGY RING 



The 4-dimension coordinate vector is X = (x T , 1 T ) T , where the horizontal coordinate vec- 
tor is x = (x, x') T and the longitudinal coordinate vector is 1 = (ct,5) t . The cr coordinate 
instead of the phase coordinate <f> = (p s — ^cr, where <p s is the synchronous phase, h is the 
harmonic number and R is the average ring radius, is used to avoid the appearance of scal- 
ing factors h/R'm the transfer matrix. Note a negative cr indicates the particle is behind 
the synchronous particle. The coordinates at the entrance and the exit of an accelerator 
component are related through its transfer matrix T which can be divided into 2x2 blocks 
M, E, F and L such that 



(1) 



In general, the coupling blocks E, F of a single component are nonzero only for dipole 
magnets. And for time-independent components (which include most common accelerator 
components except rf components such as rf cavity, rf dipoles and rf quadrupoles), the 
coupling blocks have two zero matrix elements such that 





(2) 



where e and f are 2-component column vectors and O's are zero vectors of suitable sizes. The 
zero elements in matrix E and F are consequences of the fact that the horizontal coordinates 
don't depend on the arrival time of the particles and the horizontal coordinates don't cause 
energy changes in such components. The L blocks for rf cavities and other components are 



j ri — ana lather 





(3) 



W= E^ = E — R> (4) 

where Vq is the gap voltage and E is the beam energy. The r\ parameter is related to the 
fractional phase slippage factor and is nonzero only for dipole magnets if we assume all 
particles have the same velocity c, the speed of light. 

The transfer matrix for an accelerator section is the matrix product of the transfer ma- 
trices of the sequence of components which it consists of. For any sequence of components 



not containing an rf cavity, condition Eq. (j2J) still holds. The symplecticity requirement of 
the transfer matrices of such sections is equivalent to: M and L are symplectic and 

( ° A 

MJ 2 f = e, J 2 = . 5 

\-l o) 

The e vector and 77 parameter for a section from point 1 to point 2 can be written in integral 
forms 



e 21 = / M(s 2 \s) I , (6) 




V21 = ~ / e (s 2 |s)— , (7) 



where s is the arc length along the reference orbit, s 2 |s means "from s to s 2 " , p is the 
bending radius and e 1 is the first element of e. 

The one-turn transfer matrix with an rf cavity in the ring can be derived in the following 
way. Suppose the rf cavity is located at point 2 and we want to calculate the transfer matrix 
at an arbitrary point 1. The transfer matrices between point 1, 2 and at the rf cavity are 



M21 E31 \ / M12 Ei 2 \ I I 

-21 — i I J- 12 — , J-rf 



(8) 



F21 L21 / y F 12 L12 / y L r f 

where M21 and M12 are the horizontal transfer matrices between the points, I is the 2x2 
identity matrix and L 21 and L 12 differ from the 2x2 identity matrix by their (1, 2) elements 
of 7721 and 7712, respectively. The one turn transfer matrix at point 1 is 

, Mi Ei , 

Ti 2 T rf T 2 i = , (9) 



Fi L 



(10) 



which can be expressed as 

M 12 M 2 i + Ei 2 L r fF 2 i M 12 E 2 i + Ei 2 L r fL 2 i 
Fi 2 M 2 i + Li 2 L r fF 2 i Fi 2 E 2 i + Li 2 L r fL 2 i 
We first consider the case when there is no rf cavity or the cavity is turned off. This 
corresponds to w = and L rf = I. Consequently E x and Fi satisfy conditions Eqs. (l2"f5T) . 
We intend to introduce a transformation X = UX n to block diagonalize the new transfer 
matrix T n = U _1 TU. It can be shown that this is achieved by 



where 




Dx = (I-MiPEi (12) 

and is its symplectic conjugate [7]. We may introduce a column vector di = (Di,D[) T 
so that Di = (0, di). Eq. f|T2|) defines the dispersion functions D% and D[ at point 1. One 
can also show that 

E 21 = D 2 -M 21 D! (13) 
from T 2 T 2 i = T 2 iTi. The new transfer matrix is found to be 

with Li,„ = ( 1 n I , (14) 

where Li „ is the new longitudinal transfer matrix and fj is a constant of the ring given by 

I Dds 

r) = -f—, (15) 

which is related to the usual momentum compaction factor ceo by fj = —2irRao. We have 
shown that 

f\ = rji - Hi sm2nv x , (16) 

where u x is the betatron tune, r\\ is the (1, 2) element of Lj and the 7i-function along with 
its associated phase parameter are defined by 

U = ^[D 2 + (aD + pD')% (17a) 

* = tan " (17b) 

where a, j3 are Courant-Snyder parameters. We may define the fractional phase slippage 
factor on a section between point 1 and point 2 by 

f 32 Dds 

V2i = - • (18) 

J 81 P 

It has been shown that 



V21 = V21 ~ sin(^ 2 i + xi- X2), (19) 

where ip2i = ds//3 is the betatron phase advance from point 1 to point 2. 



The new coordinates at any location after transformation are related to the original 
coordinates by 

x = Mx„ + 5 n d, (20a) 
cr = cr n + Dx' n — D'x n , (20b) 
S = S n . (20c) 

We recognize x n = (x n , x' n ) T are just the betatron coordinates. The momentum deviation 
coordinate is not changed by the transformation. But the longitudinal phase coordinate is 
different by 

Dx' n -D'x n = -V^Wcos(V'-x), (21) 
where J x is the horizontal betatron action variable and ip is the phase variable \^\ - 

III. MATRIX FORMALISM OF SYNCHROBETATRON COUPLING 

When the rf cavity is turned on, its longitudinal transfer matrix deviates from the identity 
and is now 

/ n n \ 

(22) 

Inserting it to Eq. (1101) we get the one-turn transfer matrix 

(M° E° \ 
1 1 > (23) 

where the superscript denotes quantities when the rf cavity is off and 

E12WF21 E12WL21 
L 12 WF 2 i L 12 WL 2 i 




/ Mi 


Ei 1 






Li 1 





(24) 



We may apply the procedure in Ref. |7] directly to block diagonalize matrix T in Eq. (1231) . 
However, it is easier to relate the elements in the transformation matrix to the well-known 
parameters such as dispersion functions and rf parameters in the following way. We first 
apply the transformation described in the previous section and then apply the procedure of 



Ref. [7 1 to decouple the new transfer matrix. After the first transformation, the off-diagonal 
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blocks of the transfer matrix are small because the main dispersion effect is separated. 
Hence in the second transformation we can apply some approximations to derive explicit 
expressions for the transformation matrix and the final transfer matrix. 
After the first transformation, the transfer matrix is 

T n = U- 1 TU = T°+ W f n , (25) 

where 

T o = u-i^u = M ? and f n = IT 1 * U = ( Mn E ™ I . (26) 
V Kn J \F n L 

n 

Following Sagan- Rubin [7], we want to find the transformation matrix 



V = _ (27) 




to block diagonalize the matrix T n . According to the Sagan-Rubin procedure, we define 

H = w(E n + F+) (28) 

and let 

4|I H II , x 

k = ^ ? , 29 

Tr [M n — L n ] 

where M n and L n are diagonal blocks of T n , the 7 parameter and C matrix are then 



7 = wI+l/_L C = H . (30) 

W 2 2V1 + K' 7v ^T^Tr [M n - L n ] V ; 

The elements of the H matrix have been calculated to be 

H u = -2w^f3{H 2 shi7rz/ x .cos(7rz/ x . - tp 12 - Xz), (31a) 

H12 = -w^f3iH 2 {fi sin('0i2 + X2) ~ 2^12 sin irv x cos(nv x - ip 12 - X2)), (31b) 
\h 2 

H 21 = -2wJ — S3nm/ X [aan(m/ X - ip 12 - X2) ~ «i cos(vr^. - ip 12 - X2)), (31c) 

V Pi 
\u 2 _ 

H 22 = w\ — [r/(ai sin(^i2 + X2) - cos(^i 2 + X2)) - 27712 sin 70/,; 
y Pi 

(ai cos{-kv x - ip 12 - X2) ~ sin(7rz/ x - ip 12 - Xt))]- (31d) 
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It follows that 

||H|| = w 2 ^ 2 sm2vrz/ a .. (32) 

Eqs. ()27|l29H32l constitute an analytic form of the decoupling transformation for the general 
case. In the following we simplify the expressions for the off-resonance cases in which 
Tr [M n — L n ] is not close to zero. To this end we observe that w is usually a small quantity 
(e.g., w = 0.008 m -1 for SPEAR3) and k is on the order of w 2 so that + k « 1 + 3k/8 = 
1 + 0(w 2 ). To first order of w we have 

C = 7 r, (33) 

2(cos27rz/ x — cos2ttu s ) 

where we have used Tr [M n ] = 2 cos 2irv x and Tr[L n ] = 2 cos 2irv s by definition. The parame- 
ter u x is related to the unperturbed betatron tune v x by cos 2ixv x = cos 2irv x + |u>7i 2 sin 2-kv x 
since it has been shown that Tr[M n ] = 7i 2 sin 2itv x . To guarantee the symplecticity of the 
new transfer matrix, we must have 



' X ' l|C|1 = ^ - 4(008 2^ -cos 2^)2- (34) 



The transfer matrix for the decoupled coordinates may be written as 

VEd 
L d 



. M d . 

V^V = , (35) 



where M rf and L d are given by 

M d = 7 2 M n - 7 (CF n + E„C+) + CL„C + , (36) 

L d = 7 2 L n + 7 (F n C + C + E n ) + C+M n C. (37) 

For both of the above equations, the last three terms are on the order of w 2 . To second 
order of w, we have 

M d « Ml + wM n + w 2 Mf\ (38) 

L d ^ Ll„ +wL n + w 2 L^. (39) 

Explicit expressions have been developed for M n and L n . The elements of M„ are 

Mn,n = 7Y 2 sin(^i2 + X2)(cos(^2i - X2) + «i sin(^2i - X2)), (40a) 

M n ,i2 = H2P1 sin(^i2 + X2) sin(^2i - X2), (40b) 

M n ,2\ = -z-(cos(ip2i ~ X2) + at sin(^2i - X2))(cos(t/>i 2 + X2) - «isin(^i2 + X2)),(40c) 
Pi 

Mn,22 = ^2sin(^2i - X2)(cos(^i2 + X2) - «i sin(^i2 + X2))- (40d) 



The additional terms wM„ amounts to changes to the Cour ant- Snyder parameters. For 
example, the change to (3\ to first order of w is 

= ™ U f X (cos(2 7 rz/ a: - 2^ 12 - 2 X2 ) ~ cos2nv x ). (41) 

Also, knowing the traces Tr[M^ ] = — Tr[L^ ] = —fjTi.2 cot ttu x /2, we find the total betatron 
tune change to second order of w 

Av x = -—^+ / a a . 42 

47T 167T Sin 7T^ X 

From the expressions of H and M n , it is clear that when the rf cavity is located in a 
dispersion-free region and thus Ti.2 = 0, there is no dynamic consequence from the coupling 
between horizontal and longitudinal motions. 
The matrix L n is given by 

(43) 

The unperturbed longitudinal transfer matrix L n = n + wL n describes the plain longitu- 
dinal motion. The corrections due to synchrobetatron coupling is on the next higher order 
of w. The matrix L n can be Cour ant- Snyder parametrized so that 

/ cos $ s + a s sin $ s (3 S sin $ s \ 
Ln = , (44) 

\ -7s sin $ s cos $ s - a s sin $ s J 

with /5 S 7 S = 1 + a 2 s and $ s = — 27rz/ s , where a negative sign is chosen to make f3 s always 
positive. The unperturbed synchrotron tune is given by 




The unperturbed longitudinal Courant- Snyder functions are 

1 W(f) 12 - 772l) n (f} + Wfj 12 f) 2 l) W 

a s = -77 > Ps = ^ , 7, = -r-= , 46 

2 sin 2txv s sin 2iru s sm 2irv s 

which are equivalent forms of those found in Ref . [y| . It is noted that 7 S is a positive constant 
and p s is positive but slightly varies around the ring. Without coupling, the longitudinal 
coordinates (cr n , 5 n ) are related to these parameters, the longitudinal action variable J s and 
phase variable ip s by 



2 J 

cr n = \/2J s jf s cos tp s , 5 n = - Jj^(a s cos^ s + simp s ), (47) 
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from which we obtain relations between the rms bunch length, rms momentum spread and 
the longitudinal emittance e s 

°l Tn = P s e s , °\ n = l s e s - (48) 

Note that the bunch length varies with location because particles with different momentum 
deviation experience different longitudinal phase slippage. The bunch is longest at the rf 
cavity and can be shortened by a maximum of it 2 v1/2 part of the original length at half way 
across the ring from the cavity. For fast ramping synchrotrons, it can be as large as 5% , 
assuming v s = 0.1. The synchrobetatron coupling should slightly change the synchrotron 
tune and the longitudinal Courant- Snyder functions given above. 

The decoupled coordinates X d = (xj, lJ) T are related to the original coordinates X by 

I 7 I-D 1 C+ C + 7D1 \ ( x d \ 
X = UVX d = 7 • ( 49 ) 

^_C+- 7 D+ 7I-D+C j \ U J 

The betatron coordinates (x^, 1^) T are related to the normal modes by the transformation 
matrix V. Since all four elements of the C matrix in V are nontrivial in general, the hori- 
zontal betatron coordinate x n depends on the longitudinal phase, as pointed out in Ref. 
This is a natural consequence of the synchrobetatron coupling. In fact, the longitudinal 
coordinates cr n , 5 n also depend on the horizontal betatron coordinates. In terms of the 
betatron coordinates X n , the phase space volume occupied by the beam tilts across the 
horizontal and longitudinal subspace. Particles flow in and out between the two subspaces. 
However, the total phase space volume is preserved. The reason for Ref. jfij] to suggest that 
the total phase space volume is not preserved is because it didn't fully consider the coupling 
effects on the horizontal betatron motion. 

With Eq. (j4"9l) we can decompose the 4-dimension motion to the normal modes for any 
given initial condition. It also allows us to derive the effects of synchrobetatron coupling on 
the beam sizes a x and a CT . Following the single mode analysis of Ref. [zj], we first consider 
the case when only the betatron mode, or mode a, is excited so that 

x 
x 1 

where (x a , x' a ) are the betatron normal mode coordinates given by 



/ X a ' 


!• 1 








T 



(Tl-DxC*) , =-(C + + 7 Dr) , (50) 




2 J 

x a = ^2J a f3 a cosifj a , x' a = - J j^{a a cosip a + sinipa). (51) 
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It follows that 



b^Wa + D(C 21 ^/K+ Cn-%)} cos + sin (52a) 



< T -WaD' + a a D) + {(3 a C 22 + a a C l2 )}^=co^ a - { ° 12 tJ D) sim/v (52b) 



The bunch width and length can be derived from the above equations by integrating over 
the bunch distribution. Here we consider only the off-resonance cases and assume cos v s w 1 
so that C ~ H/4sin 2 7ri/ x . We get 



= PaCa-e a — : sm^^ - ^12 - X2) + ea— — -2 {D - -Pir}C0tiri> x ),{b3) 

sin 7ru x 4 sin 7rz/ x 2 

where e a is the emittance of the horizontal betatron normal mode. We see that the term on 
the order of w varies around the ring. The leading non-varying correction term is on the 
order of w 2 . Since the 0(w 2 ) terms are very small, we will drop them in the following for 
brevity. Similarly for the bunch length we obtain 



a 2 CTa = e a U\ - e a W , Hl ^ 2 [rj cos(^i 2 + X2 ~ Xi) ~ 2r] 12 sm-KU x sin(7rp x - ip 12 -X2 + X1)], 
2 sin 7ru x 

(54) 



where is the 7i-function at point 1 evaluated with the perturbed Courant-Snyder func- 
tions. 

The same analysis can be applied to the synchrotron mode, or mode b, to obtain its 
contributions to bunch width and length, which are given by 



2 n 2 , wr)D^/PJi 2 

a X b = e blbD +e b j b . cosTTU x sm('Ki' x -ip 12 -X2), (55) 

2 sin 7rz/ x 



a 



CT,b 



tbPb - 6b7&— cos(-kv x - ip 12 -X2 + Xi)- (56) 



sm it u x 

Note that 7^ is an 0(w) factor. The momentum spread is given by 



2 w 2 H 2 2 w 2 r]H 2 sin 2irv x . . 

a S,a = e «7-l > a S,b = e blb ~ Cblb — . 4 • (57) 



X 



So to first order of w we have a 2 = e^b- 

The contributions to cr 2 from mode a and 6 simply add up i.e. 



^ = < + cr 2 , (58) 
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which applies to all four coordinates. Eqs. (I53H57P show that because of synchrobetatron 
coupling, the longitudinal motion affects the transverse beam size and and the transverse 
motion affects the longitudinal beam size. The leading correction term is usually a small 
term on the order of w and varies with the horizontal betatronphase advance around the 
ring. It is noted that we have recovered Shoji's result in Ref. [5( by the first term in Eq. 
(|5^|) which is independent of w and indicates that the bunch length varies from location to 
location according to the local 7i-function. This is a consequence of the uneven distribution 
of the path length effect. 



IV. CLOSED-ORBIT CHANGE DUE TO ENERGY GAIN AND LOSS 

The finite energy gain at the rf cavity and the radiation energy loss around the ring are 
kicks to the momentum deviation coordinate. These kicks are transfered downstream and 
affect all other coordinates. Therefore the closed orbit of the beam is changed. To study this 
effect in the matrix formalism, we extend the coordinate vector to X e = (x, x' , cr, 5, 1) T sj. 
The corresponding transfer matrix is then 5x5. The fifth element is included to describe 
the kicks received by the particle, namely the equation 

X 2 = T 21 Xx + g 21 (59) 

will now be written as 

X e 2 = T^X?, T e 21 = ( T q 21 ^ j (60) 

where g 2 i is a 4-vector which represents the kick-induced shifts of phase space coordinates 
across the accelerator section from point 1 to point 2. The closed orbit X c is given by the 
fixed point X^; = (X^T, 1) T of the extended one-turn transfer matrix T e sj], i.e., T e X^ = XJ;, 
which yields 

X c =(I-T)- 1 g, (61) 

where the g vector contains the first four elements of the fifth column of T e and represents the 
coordinate shifts after one turn when a particle starts from a point with initial coordinates 
of all zeros. We will derive analytic forms for g and (I — T) _1 below. 
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For an rf cavity, there is an energy kick 

g rf = (0,0,0, e) T , 



AE 

IT' 



(62) 



where e denotes the sudden change of momentum deviation at the rf cavity. For a dipole 
magnet, the energy kick due to radiation energy loss will propagate to the other coordinates 
and cause finite changes to them. For example, the pure sector dipole with bending radius 
p and bending angle 0l has 

/ p(9 L -sm9 L ) \ 
1 — cos 9 L 



Sdipole 



2tx ! 



p(l-caaO L -\0 2 L ) 
6 L 



(63) 



/ 



where —69l/2ti is the momentum deviation change on this dipole magnet. Here we have 
assumed e is the same as the momentum deviation gained at the rf cavity. It is seen that 
the changes of x, x' and cr are on the order of Q\, 9\ and 9\, respectively. Hence when Q L 
is small, we may neglect these changes for a single dipole magnet. The g vectors are zeros 
for other accelerator components which don't cause energy gains or losses. 

The one turn extended transfer matrix can be computed as usual. For an arbitrary point 
1 , we have = T^T^T^, from which we obtain 



gi = (Ti 2 g 21 + gi 2 ) + T 12grf + T 12 (T rf - I)g 21 . 
Simple calculations show that 



Tl 2 g rf 



(64) 








( di - M 12 d 2 ^ 












= e 






1 ('") ) 


W 







(65) 



The T 12 g 21 + g 12 term has nothing to do with the rf cavity. So it does not depend on the 
location of point 2. In fact it represents the changes of coordinates in one turn for a particle 
with initial coordinate (0, 0, 0, 0) T at point 1 when the rf cavity is turned off. Since energy 
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losses occur in dipoles, we get a summation over all dipoles in the ring 

/ \ 

Tl 2 g21 + gl2 = 2^ 



=1 \ Fij Lii 



/ 



o 

Noting that — di — M^dj and, assuming all bending radius are equal so that e« 
eAsj/27rp, the summations can be turned to integrals to obtain 



(66) 



/ 



Tl 2 g21 + g 
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d 1 -&M{8 1 \8)d{8)g 



ds 
2-rrp 



\ 



\ 



J 



(67) 



where we have used Ei e « = e an d the integral § x starts from point 1, completes one rev- 
olution and ends at point 1. We may simplify the expressions by defining new functions 



S{s) 
C(s) 



s+C 



\jH(s') sin(Vw + Xs', 



s+C 



\/H{s') cos(V> ss' + 



ds' 
ds' 



Xs 



/C(s) = S 2 (s) +C 2 (s), f(s)=tan 



c( s y 



(68a) 
(68b) 
(68c) 



where ip ss > is the betatron phase advance from point s' to point s and C is the ring circum- 
ference. Then we can write down 



ds 



%. 

Ctl 1 



Si 



and 



M(si|s)d(s. , . , 

r \ ,/~R7 ./~R7 I \ 1 



, , s ds 
V(si\s) — 
i 2vrp 



(69) 



^(si|s)— + y/HiKi sin(£i - xi). 



It can be shown that 



_ ds 1 _ 
1 ^ 1|S) 2^ = 2^ 



(70) 



(71) 



Since the first three elements of the g vector have terms on the order of 0(e) and the 
contributions of the third term of Eq. (|64p to these elements are on the order of 0(we), we 
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will simply drop these contributions. However, the fourth element of this term is the leading 
term and is not negligible. It is easy to show that 

ds 



rs2 

9\ = wgti = / r](s 2 \s) 

J Si 



with the definition of 



2vrp 

/S2 ^ s 
Vi s 2\s)- we^jH 2 K, 2 i sin(^2i - X2) (72) 
.1 2lT P 

I-S2 

S21 = I a/H(s') sin(^ S2S ' + Xs')o — ' ( 73a ) 

fS2 

C21 = / ^fH{s>) COSTS' +Xs>) 7T- > ( 73b ) 



'si 



£21 = S 2 21 + C 2 2 l5 61 = tan" 1 (73c) 

If we write down the vector gi in a form g x = (g^, gf ) T with 2-component vector g x and g^ 
representing the horizontal and longitudinal displacement, respectively, we have 

./«= I • • ; 11 + ^ (74a) 

Ci - V^2Cos(V'i 2 + X2) 




9i/e = 7712 - -77 + \fHiH2 sin(^i2 + X2 - Xi) - \/^i^i sin (6 - Xi), (74b) 



—2 

9i/z = ~ w ip~ ~ wy/Wn sinfei - X2), (74c) 

where gl, gf are the two elements of g; = {gi,gf) T and in obtaining the first term of Eq. 
(I74c| we have assumed the dipoles are distributed around the ring uniformly. 

Radiation damping changes the transfer matrix T by a correction term on the order of e 
so that it is no more strictly symplectic. However, since e is usually small and the correction 
to the closed orbit due to this effect is on the order of e 2 , we will not consider the radiation 
damping effect. To work out the matrix inversion for (T — I) -1 , we will make use of 

(T - I)" 1 = U(T n - I)- 1 !!" 1 , (75) 

since the off-diagonal blocks of T„ are on the order of 0(w). With matrix T n as found in 
Eq. (I25I) . we have shown that 

(TWr 1 = I S (7G) 
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with 



a = (M°-I)~\ (77a) 

b = -(M5-I)- 1 | ° E » | , (77b) 

E n 

(M? - I)- 1 , (77c) 



P 21 F 22 

■ n n 



d = (L^-I)- 1 , (77d) 
where we have dropped all terms on the order of 0(w) or higher and we have 

(L.-D- 1 = -I"' 2 ' - +Ml2 V (78) 

' V 1 - ,iu I 

We then proceed to obtain the blocks of (T — I) -1 using Eq. (1751) and finally we get the 
closed orbit with Eq. (16T1) . The results are given by 

cos(7r^ - - X2) + 7T-- cos{tcv x - £1) +eD 1 {- — ), (79a) 



x 



c 



2siD.Tvu x ' 2sin7rz/ x . 2 77 



x 



^ (ai cos (71-1/3 - ^12 - X2) - sin(7rv x - ^12 - X2)) 



2 sin 7r^ y (3 



6 / X ~~ 

7r~. a -^■(a 1 coa(m/ x - £1) - sin^z^ - &)) + eDj_(- - — ), (79b) 

2 sm 7^ y pi 7 2 r; 

cr c = 6 21 } 2 + — [y/HiHz cos(irv x - ifa - \2 + Xi) + cos(tt^. - ^12 - X2 

277 2 sm 7rz/ x L 



+£1) - v 7 ^1^1 cos(7r^ + xi ~ 6) - ^2 cos7rz/ x . + e \/Ji^2\ sin(£ 2 i - X2), (79c) 

5 c = 4-%. (79d) 
2 77 

Note again that we have dropped all terms on the order of w or higher so that the results 
are valid only for off-resonance cases. Obviously we have recovered Eqs. (23-24) of Ref. 6| 
as the first term in Eq. (I79aj) . The second term comes from the propagation of radiation 
energy losses. The third term comes from the energy variation around the ring. It is worth 
noting that the cr c orbit is zero at the rf cavity. 

So far we have considered only one rf cavity in the ring. However, the same analysis 
can be easily applied to more cavities. In fact, if we neglect the interaction between the rf 
cavities, which corresponds to higher order terms of the Wi parameters, the one-turn transfer 
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matrix is 



Tx = T° + J^TWiTa, (80) 

i 

where the summation is over all cavities. Eq. (180]) means the total effect of all cavities is 
the linear superposition of the single cavity effects. It is then straightforward to modify the 
results of the one-cavity case for multiple-cavity cases. 



V. SIMULATION 



In this section we will verify the theory developed in the previous sections by comparing 
the results to simulations with the accelerator modeling code AT 12j. We use the machine 
model of the SPEAR Booster for the calculation. The SPEAR Booster is chosen because 
it has appreciable dispersion functions at the rf cavity. The model consists of 20 periods of 
FODO lattice over a circumference of 2ttR = 133.8 m. The bending radius is p = 11.82 m 
for all dipoles. The extraction energy is E = 3 GeV. The rf frequency is f r { = 358.533 MHz 
and the harmonic number is h = 160. At extraction, the rf gap voltage is V T { = 0.8 MV and 
the one-turn radiation energy loss is Uq = 0.60 MeV. In the simulation, we will consider it 
as a storage ring running at its extraction energy. The betatron tunes are u x = 6.16 and 



4.23 for the model. At the rf cavity, the Courant-Snyder functions are a<i = —0.72 and 



02 = 2.09 m and the dispersion functions are D 2 = 0.27 m and D' 2 = 0.06. The horizontal 
chromaticity is C x = —7.9. 

We first present a numerical example to verify the transformation Eq. (1491) . The injection 
point of the SPEAR Booster is located in the ring opposite to the rf cavity where we have 
q;i = 0.02, f3\ = 1.54 m, D\ = 0.24 m and D[ = 0.01 and the betatron phase advance from 
the rf cavity to this observation point is ipi2 = 0.16 rad modulo 2n. The synchronous phase 
is set to <p s = tt and hence w = 0.0020 m -1 . The one-turn transfer matrix for betatron 
coordinates X„ at this location is 



/ 0.550532 1.306266 -0.000447 0.001018 \ 
-0.547728 0.516813 0.000129 -0.000294 

0.000396 -0.000910 0.995433 -4.548267 

\ -0.000174 0.000399 0.002004 0.995433 / 



(81) 
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(82) 



The off-diagonal blocks represents the synchrobetatron coupling effects. We then apply the 
second transformation V as determined by Eqs. (1271 I3TI I33H341) to obtain the final transfer 
matrix 

/ 0.550533 1.306265 -4.6 x lO" 10 1.05 x 10~ 9 \ 
-0.547728 0.516813 1.3 x lO" 10 -3.0 x lO" 10 
4.1 x 10~ 10 -9.4 x 10^ 10 0.995432 -4.548266 
\ -1.8 x lO" 10 4.1 x 10~ 10 0.002004 0.995432 / 

Elements of the off-diagonal blocks are reduced to the order of magnitude of 10 -10 from the 
original value of 10~ 4 , indicating the high precision of the approximations we have made. 

The betatron detuning due to synchrobetatron coupling (Eq. (T4"2l ) is checked against AT 
by comparing the betatron tune from AT tracking to the analytic calculation with various 
rf voltage. The initial condition for AT tracking is (1 x 10~ 5 , 0,0,0, 1 x 10~ 6 ,0) T . Small 
deviations of x and 5 are chosen to avoid significant nonlinear detuning effects. The result 
is shown in Figure [TJ 
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FIG. 1: The synchrobetatron coupling induced betatron detuning Az^ x obtained by AT tracking 
is compared to analytic calculation with Eq. (|42p . The rf gap voltage is varied from 0.2 MV to 
3.2 MV. 



If radiation is turned on in simulation, the finite energy gain at the rf cavity will cause 
changes to the closed orbit as studied in the previous section. To calculate the closed orbit 
change, we need to calculate the functions S, C, JC and £ defined in Eq. fl68l) . It is found 
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that numerical integration through a dipole does not make much difference from simply 
transporting the coordinate shifts at the dipole exit to the observation point. In Figure 
[2] we show functions S and C for the SPEAR Booster ring calculated with both methods. 
The one-turn coordinate shifts, or the g vector is plotted in Figure [3J There are two curves 




s (m) 

FIG. 2: (Color online) Function S and C as defined in Eq. (|68|) for the SPEAR Booster at 3 
GeV. The curves labeled "Int" are from numerical integration with the Trapezoidal rule by cutting 
each dipole into 100 slices. The "Non" curves are obtained by calculating the coordinate shifts at 
the exit of each dipole, transporting to the observation point and summing up contributions of all 
dipoles. 

in each plot, one is from Eq. (|74|) . the other is from direct matrix multiplication with the 
transfer matrix of each accelerator element given by AT and the coordinate shifts of each 
element obtained by tracking through it with zero initial coordinates in AT. 

We then compare the closed orbit changes obtained with three different ways: (1) AT 
(using the function findorbit6), (2) direct matrix calculation with Eq. (I6T1) and (3) the 
analytic solution in Eqs. ()79aH79dl) . The results are shown in Figure HI Good agreement 
between the three curves are seen, verifying the analysis in the last section. The S and 
C functions are calculated with simple summations in the analytic approach. It is worth 
pointing out that all the three terms in Eq. (179aj) are important in determining the closed 
orbit change for x c . 

The same closed orbit calculation is done for the SPEAR3 ring. We only show the closed 
orbit changes in Figure Note dispersion at the rf cavity for the SPEAR3 ring is zero. So 
the terms in Eqs. ()79aH79d|) involving 7i 2 have no contribution. 
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FIG. 3: (Color online) The one-turn coordinate shifts due to radiation energy losses and finite 
energy gains at the rf cavity for the SPEAR Booster at 3 GeV. The four plots are x, x', ct and 5 
shifts, respectively. 




s (m) s (m) 

FIG. 4: (Color online) The SBC induced closed orbit changes calculated in three different ways for 
the SPEAR Booster at 3 GeV. The changes to x c , x' c , ct c and 5 C are shown in the four plots. The 
three curves are from AT ("AT"), matrix inversion ("Mat") and the analytic solution ("Ana"), 
respectively. 
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VI. CONCLUSIONS 



In this study we fully analyzed the linear synchrobetatron coupling by block diagonalizing 
the 4x4 horizontal-longitudinal transfer matrix. We found the transformation between the 
usual (x, x', ct, 5) coordinates and the normal modes and the transfer matrix for the normal 
modes in analytic forms in terms of the Courant- Snyder functions, dispersion functions and 
the rf voltage slope parameter. This enabled us to predict the 4-dimensional motion of 
a particle knowing only the initial condition and those common parameters. The effects 
of synchrobetatron coupling on the horizontal betatron motion, including changes to the 
Courant- Snyder functions and the betatron tune are also presented. We then studied the 
beam width and the bunch length under synchrobetatron coupling. We readily recovered 
Shoji's result of dispersion- dependent bunch lengthening [5|. We found that the beam width 
and the bunch length slightly oscillate around the ring with the betatron phase advance 
measured from the rf cavity. We pointed out that the bunch length varies around the ring 
due to phase slippage, a fact that is often overlooked. We also pointed out that the phase 
space volume is preserved under SBC when not considering radiation induced diffusion and 
damping. 
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FIG. 5: (Color online) The SBC induced closed orbit changes calculated in three different ways 
for the SPEAR3 ring. The changes to x c , x' c , ct c and 5 C are shown. The three curves are from AT 
("AT"), matrix inversion ("Mat") and the analytic solution ("Ana"), respectively. 
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Following Ref. Q], which studied the horizontal closed orbit changes due to finite energy 
gains at rf cavities, we fully explored the problem by looking for a closed orbit in the 4- 
dimension phase space, considering both finite energy gains at rf cavities and radiation 
energy losses in dipole magnets. We recovered the horizontal closed orbit change result in 
Ref. p and found additional terms due to energy losses. 

We carried out simulations with the accelerator modeling code AT 12j to verify the 
analysis. We found that the block diagonalization transformation matrix had high precision 
despite the approximations we made to get the analytic solution. The closed orbit results 
of AT also agreed well with the theory. 
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